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ABSTRACT 

We analyze the data of the gravitational microlensing survey carried out by 
by the MOA group during 2000 towards the Galactic Bulge (GB). Our obser- 
vations are designed to detect efficiently high magnification events with faint 
source stars and short timescale events, by increasing the the sampling rate up 
to ~ 6 times per night and using Difference Image Analysis (DIA). We detect 
28 microlensing candidates in 12 GB fields corresponding to 16 deg^. We use 
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Monte Carlo simulations to estimate our microlensing event detection efficiency, 
where we construct the /-band extinction map of our GB fields in order to find 
dereddened magnitudes. We find a systematic bias and large uncertainty in the 
measured value of the timescale tEout in our simulations. They are associated 
with blending and unresolved sources, and are allowed for in our measurements. 
We compute an optical depth r = 2.59lQ;g4 ^ towards the GB for events 
with timescales 0.3 < < 200 days. We consider disk-disk lensing, and obtain 
an optical depth Tbuige — 3.361q;8} x 10~^[0. 77/(1 — /disk)] for the bulge compo- 
nent assuming a 23% stellar contribution from disk stars. These observed optical 
depths are consistent with previous measurements by the MACHO and OGLE 
groups, and still higher than those predicted by existing Galactic models. We 
present the timescale distribution of the observed events, and find there are no 
significant short events of a few days, in spite of our high detection efficiency for 
short timescale events down to ~ 0.3 days. We find that half of all our de- 
tected events have high magnification (> 10). These events are useful for studies 
of extra-solar planets. 

Subject headings: dark matter — Galaxy:halo — gravitational lensing 



1. Introduction 

Following the suggestion of Paczyhski (1991) and Griest et al. (1991), several groups have 
carried out microlensing surveys towards the Galactic Bulge (GB), as seen in Baade's window. 
It is now well understood that these observations are useful for studying the structure, 
dynamics and kinematics of the Galaxy and the stellar mass function as the event rate and 
timescale distributions are related to the masses and velocities of lens objects. 

The amplification of a microlensing event is described by (Paczyhski 1986) 

= r^— , (1) 

where u is the projected separation of the source and lens in units of the Einstein radius Re 
which is given by 

i?E(M, x) = ^^D,xil-x), (2) 
where M is the lens mass, x = D\/ Dg is the normalized lens distance and D\ and Dg are the 
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observer-lens and the observer-source star distances. The time variation oi u — u{t) is 



where P, to, t-E = Re/vi and ft are the minimum impact parameter in units of Re, the 
time of maximum magnification, the event time scale and the transverse velocity of the lens 
relative to the line of sight towards the source star, respectively. Prom a light curve, one can 
determine the values of /3, to and ^e, but not the values oi M, x or vt- 

Our MOA (Microlensing Observations in Astrophysics) group started observations to- 
wards the GB in 1999. Prom 2000 we introduced the Difference Image Analysis (DIA) 
(Crotts 1992; PhiUips & Davis 1995; Alard & Lupton 1998; Alard 2000; Alcock et al. 1999, 
2000; Wozniak 2000; Bond et al. 2001) which is able to perform better photometry than the 
traditional DoPHOT (Schechter, Mateo & Saha 1993) type analysis in crowded fields at any 
place even where no star was identified. 

To date, hundreds of microlensing events have been detected towards the GB by the 
OGLE (Udalski et al. 1994, 2000; Wozniak et al. 2001) and MACHO collaborations (Alcock 

et al. 1997a, 2000). They estimate the microlensing optical depth towards the GB to be 
3.3 ± 1.2 X IQ-s from 9 events by DoPHOT analysis, 3.9+}:2 x IQ-^ from 13 events in a clump 
giant subsample from DoPHOT, and 3.23lg5Q x 10~^ from 99 events by DIA respectively. 
Popowski et al. (2000) and Popowski (2002) estimate values of 2.0toA x 10"® and 2.23l[]:| X 
10~^ respectively from MACHO data. These values are all more than twice those expected 
from existing Galactic models, which are somewhere around 0.5 ~ 1.0 x 10~^ (Paczyhski 
1991; Griest et al. 1991; Kiraga & Paczyhski 1994). This suggests that the standard models 
of the Galaxy need to be revised. To explain the high optical depth, a number of authors 
have suggested the presence of a bar oriented along our line of sight to the GB (Paczyhski et 
al. 1994; Zhao, Spcrgcl & Rich 1995), and have adopted various values of the bar orientation 
and mass (Paczyhski et al. 1994; Peale 1998; Zhao & Mao 1996). Microlensing observations 
towards the GB therefore appear useful for characterizing the mass and inclination of the 
bar. 

Popowski et al. (2000) and Popowski (2002) raised the possibility of a systematic 
bias in the optical depth due to the difficulties of measuring tE associated with blending 
and unresolved sources. When the actual source base-line flux is unknown, Ie and j3 arc 
degenerate in relatively low signal-to-noise ratio (S/N) events (c.f. Han 1999; Bond et al. 
2001;Gould & An 2002). The optical depth may be estimated by using red clump giant 
stars to avoid the bias, or by other methods (e.g. Gondolo 1999; Kerins et al. 2001). In 
this paper, we quantify the bias by using Monte Carlo simulations and take itinto account 
to estimate the optical depth. 
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Our observations are designed to detect efficiently high magnification events with faint 
source stars for the study of extra-solar planets and surface-transit events, by increasing 
the sampling rate up to ~ 6 times per night (note that the sampling rate of other projects 
is typically once per night). Our observations are consequently fairly sensitive to short 
duration events, i.e., events caused by smaller mass lenses. This could lead to a different 
optical depth estimate from previous studies if there is a significant contribution of low mass 
objects such as brown dwarfs to the microlensing optical depth. Thus the MOA observations 
can constrain the contribution of the low-mass population and also the structure of the GB. 
For this purpose here we present the results of MOA observations. 

In this paper, the results of the DIA analysis of data towards the GB taken by MOA 
in 2000 with DIA are presented. The analysis is aimed at finding how efficiently we can 
detect high magnification events and short timescale events, and at estimating the optical 
depth towards the GB. In § 2 we describe our observations. § 3 is devoted to the analysis 
method. In § 4 we describe the event selection process and results. In § 5 wc make an J-band 
extinction map of our GB fields caused by dust. This is useful in estimating extinction free 
source magnitudes in the simulation. In § 6 we describe the simulation used to estimate our 
detection efficiency and the resultant microlensing optical depth. Discussion and conclusions 
are given in § 8. 



2. Observations 

We observed from the Mt.John University Observatory (MJUO) in New Zealand at 
170°27.9'E, 43°59.2'S, 1030 m ahitude. Typical sky background values are 21.9, 22.6, 21.5, 20.9 
and 19.1 magarcsec^^ in the U, B, V, Rc and Ic passbands respectively at an air mass of 
1.0 (Gilmore, private communication, 1994). Spectroscopic hours arc ~ 45%. 

We used a 61-cm telescope equipped with the large mosaic CCD camera M0A-cam2 
(Yanagisawa 2000). This has three 2k x 4k pixels thinned CCDs. The combined field of view 
of the camera and telescope is 0.92° x 1.39° (0.81 arcsecpixel^^). We used two non-standard 
wide passband Red (630 — 1100 nm) and Blue (400 — 630 nm) filters. The global seeing is 
typically 1.9 - 3.5". 

We observed 14 GB fields (18 deg^) 3 ~ 6 times per night during the southern winter 
season (from April to November) , where each field consists of three subfields corresponding 
to each CCD. The main aim of the observations was to detect high magnification events in 
order to find extrasolar planets (Bond et al. 2001). The detection probabihty of extrasolar 
planets in high magnification events is high (Griest & Safizadeh 1998). 
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Microlensing events in which the source star is fainter than the observational hmiting 
magnitude, are generally known as "Pixel Lensing" events (Gould 1996) or "EAGLE events" 
(Nakamura & Nishi 1998; Sumi & Honma 2000). These usually have high magnification. 
However, the period in which they are visible is very short, because they are visible only at 
peak of magnification. To detect this kind of events we sample each field 5 ~ 6 times per day 
mainly through the Red filter rather than both colors. Due to this high sampling rate, the 
detection efficiency for short timescale events also becomes higher than the former studies. 
A unique feature of our observations is the ability to detect small lens objects towards the 
GB. 



3. DATA Analysis 

We have used the DIA technique to provide sensitivity to Pixel Lensing or EAGLE 
events. These events cannot be detected by the conventional (fixed position) DoPHOT-type 
photometry. 

We developed our own implementation of the DIA method of Alard & Lupton (1998) 
and Alard (2000) which directly model the kernel in real space with spatial variations of the 
kernel across the CCD. This method is suitable even for crowded and poor S/N images. In 
this technique, one first registers some star positions, and geometrically aligns each "current" 
observation image to a preselected "reference" image which is recorded in good seeing, high 
S/N and low airmass. Next, the convolution kernel which is required to map the reference 
image to the current image, is calculated by using the current and reference images. The 
reference images are convolved to match the seeing and scaled to align with the current 
images. The resultant images are then subtracted and variable objects are detected. 

The images were analyzed in real-time at MJUO to issue alerts for the events (Bond et 
al. 2001). In this analysis we use the same sample of subtracted images as in the real-time 
analysis. 

We have found the effect of differential refraction, which causes a significant residual 
flux in the subtracted image (Tomany & Crotts 1996; Alcock ct al. 1999), is negligible in 
our images with the red filter when the airmass is lower than 2.0. We have not corrected for 
this effect because our data are mainly taken through the red filter. 

We have made a catalogue of all stars in the reference images by using DoPHOT. All 
saturated stars and bad pixels on the reference images are masked out the same as in the 
subtraction process. These are useful in the following analyses. 
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3.1. Identification of variable objects 

On a subtracted image, variable objects can be seen as positive or negative profiles 
depending on whether the fiux had increased or decreased relative to the reference image. 

To detect these objects, we use our implementation of the algorithm in the IRAF task 
DAOFIND, in which both positive and negative profiles could be detected simultaneously. 
There are also spurious profiles not associated with stellar variability, for example, cosmic 
rays, satellite tracks and electrons leaked from bright saturated star images. To avoid de- 
tecting these spurious objects, we apply several criteria about the statistics and profiles in 
the Analog Digital Unit (ADU) value of pixels around the peak, e.g., S/N is required to be 
larger than 3. 

Candidate objects which pass these criteria are checked against those obtained in pre- 
vious reductions of the field. If no object is cross-referenced, these new objects are added 
to the initial list of variable object positions. If an object has previously been detected, 
this is identified as the same object and the number of detections for this object is incre- 
mented. The number of detections for each object is used in the event selection cut "cutl" 
(see §4.1.1). 

3.2. Photometry 

To the list of variable objects made in the previous section, we have applied the first 
simple event selection "cutl" (see §4.1.1) to cull the number of spurious objects in this fist. 

For the objects which passed the "cutl" test, we have performed PSF profile fitting 
photometry, where the high S/N empirical PSF images (23 x 23 pixels) is made for each of 
the 32 (500 x 500 pixels) sub-regions in the reference image. Then these PSF images are 
convolved by the same kernel function used in the subtraction process, to match this PSF 
to that in each local position of each time series of frames. And the total fiux of this PSF is 
normalized to 1. 

The PSF photometry comes down to a two-parameter fit for the amplitude Upgf and base 
fine 6, i.e. apsfPSFi + h, here PSFi is the value of empirical PSF image in pixel i centered 
on the variable object. In principal, in subtracted images, the background is zero by virtue 
of the least-square process. However, in fact, there are residuals especially around the very 
bright saturated stars, bright variable stars and in poor S/N images. In our images, the 
fraction of the area where these effects occur is small but not negligible because of the wide 
wing of the PSF (our median seeing is about 2". 3). When we compared the 1-parameter and 
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2-parameter fits, we found that several faint stars around a bright variable star show similar 
variability as the bright variable star in the 1-parameter fit, but not in the 2-parameter fit. 
Fitting a 2nd "base line flux" parameter on the subtracted images provides a good zero-point 
check. So we make the base line b a fitting parameter. In the subtracted images the fiux 
has been rescaled to match that in the reference image by dividing with the scale factor oq 
which is the first coefficient of the kernel function of the frame. By this process, all the time 
series of subtracted images are photometrically calibrated. Then no extinction correction 
is required. In the fitting, the noise of each pixel (7j of the subtracted image is given from 
the actual flux fi in the current image before subtraction, taking account of the gain G in 
ADU/e~ and the scale factor of each frame. 

This photometry has been performed on all objects in the variable list in all the time 
series of images and the differential flux AF — apgf light curves stored in the database with 
their corresponding error and the square root of the reduced the PSF fltting (we refer 
to this as SdeVpsf) . 

3.2.1. Noise properties 

The error in the PSF photometry apsj obtained by fitting would be optimistic, as 
it includes only the photon noise component. So we test the properties of the photometry 
using a sample of constant stars. We have randomly sampled ~ 1000 constant stars for each 
image from our star catalogue, in which we have rejected the variable objects detected in 
§3.1. 

Noise properties are derived from the residuals of individual measurements around the 
mean flux of each object. So we use ~ 100, 000 measurements for each chip. Each residual is 
normalized by the error from the photon noise apsf for the corresponding photometric point. 
Then stars are grouped according to their brightness on the reference image, the standard 
deviation normalized by apgf in the PSF fitting, i.e. SdeVpgf. We think that the SdeVpgf of 
each measurement should be a good indicator for the systematic noise which comes from the 
non-photon noise. All residuals coming from the light curves of stars in a given group are 
merged into one distribution. 

For each group we calculate the half width of the region containing 68.3% of residuals 
'^resid-i ^ robust estimator of the width a of the Gaussian distribution. In Figure 1, we plot 
the estimated aresid as a function of SdeVpsf and these are fitted with a 3-degree polynomial. 
Here bright stars {R < —11 in DoPHOT output red magnitude, which corresponds to / < 15; 
see § 3.3) are rejected because such bright stars have other systematic deviations as discussed 
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below. The effect of such bright stars is small since, in total, they represent less than 0.4% 
of all possible source stars (/ < 23), and the pixels where such stars dominate on the CCD 
chip is less than 1%. We did not use those measurements with SdeVpsf > 10 as they are 
unreliable. The fraction of such measurements is ~ 7% out of all measurements. In this 
figure, a clear trend for these values can be seen, and the fitted curve can be used to rescale 
the error bars in order to improve the consistency of the light curves. We also show the same 
plot normalized with this fitted function in Figure 1. 

To check whether this normalization is appropriate or not, we show the distribution 
histograms of these residuals (upper panel) and normalized residuals (lower panel) in Figure 
2. As shown in Figure 2, the distribution of normalized residuals is consistent with the 
normal distribution. So our normalization seems to be reasonable. This exercise have been 
carried out for each of the CCD chips which showed that chips 2 and 3 are very similar but 
differ slightly from chip 1. So we have estimated this function for each chip. This trend 
is very similar over all fields measured with the same CCD chip. We have applied this 
normalization to all our measurements. 

We also show the Oresid of residuals (filled circle) and normahzed residuals (cross) as a 
function of brightness of constant stars in 7-band magnitude (see § 3.3) in Figure 3. As shown 

in Figure 3, as we move towards brighter stars the (Jresid increases due to systematic effects 
related to the uncertainty in the seeing and the PSF, first slowly and then faster. Alard 
& Lupton (1998) provide a possible explanation in terms of the atmospheric turbulence. 
Another possible source of this excess is the existence of some small amplitude variable stars 
in these bright stars. No correction for this effect is made in this analysis. In Figure 3, we 
find that the noise (standard deviation) in this analysis is about 40% above the photon noise 
limit. 



3.3. Calibration of Fluxes 

To obtain the source magnitude from the baseline flux /o given by light curve fitting, we 
need the transformation relation between the flux in the reference image and the apparent 
magnitude in standard passbands. We use the UBVI photometry catalogue of selected stars 
in Baade's window provided by OGLE (Paczyhski et al. 1999) that are contained in three 
MOA GB subfields (ngb2-l, ngb2-2, ngb3-3, where ngb2-l means chip 1 in field ngb2). In 
these regions, the extinction is relatively low and uniform. Our star catalogue have been 
made by applying DoPHOT to half of all the frames for each field which contain ~ 20 Blue 
images taken during the 2000 season. Then we have made the R{MOA) and B{MOA) 
catalogue of all stars in our GB fields by taking the median for each star so as to increase the 
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accuracy of the photometry and to avoid the daily difference of extinction by the atmosphere. 
Then the DoPHOT photometry i?re/ of the reference images is compared with the catalogue 
value Rmed- The differences are Rref — Rmed = —0.1 ~ +0.03 mag and spatial variations in 
one chip arc about ~ 0.02 mag. So we have neglected local differences in the chip and stored 
one offset value for each chip. These offsets are used in the following calibrations. 

Next, the R and B measurements in our catalogue of 3 MOA-OGLE overlap fields are 
compared with the V and / photometry of the corresponding star in the OGLE catalogue. 
Fundamentally the transformation to the standard V and / from MOA non-standard R and 
B depends on the color index B — R. The transformations obtained to the standard / and 
V for each chip are shown in Figure 4 and 5, and are given by following equations. 



7 = - 0.0969(5 - i?) + 26.2840 (chipl), (4) 

I^R- 0.0969(5 - i?) + 26.3331 (chip2) , (5) 

7 = i? - 0.0969(5 - i?) + 26.5937 (chip3), (6) 

V^B- 0.160(5 - 7?) + 26.350 (chipl), (7) 

V^B- 0.160(5 - 7?) + 26.800 (chip2), (8) 

y = 5 - 0.160(5 - 7?) + 26.192 (chip3). (9) 



Here the slopes were estimated by the using full data of all 3 chips, and the offsets were 
estimated individually for each chip. 

We apply these transformations to our star catalogue. As shown in Figure 4, the color 
dependence in transformation to 7 from 7? is weak, and we find that the transformation 
without a color term seems to work properly, as shown in Figure 6. The transformation 
functions can be written as follows 

7 = 7^ + 26.0923 (chipl), (10) 
7 = 7^ + 26.2210 (chip2), (11) 
7 = 7? + 26.4963 (chip3). (12) 



These equations are very useful because our data were usually taken only in Red. 
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4. Results 

All data consist of 7,200 images (~ 100 GByte) which correspond to 118 ~ 224 frames 
for each field depending on their priority. Of these, 6, 600 images (~ 92%) have been suc- 
cessfully reduced. The reduction failures occur as a result of poor data conditions, i.e., bad 
pointing, poor seeing and very low transmission due to clouds. 

In these subtracted images we find 2, 000 ~ 10, 000 variable objects in each field. And 
200 ~ 1,000 of these variable objects pass the "cutl", for which the PSF photometry is 
performed to make the light curves. For these fight curves, we apply an additional selection, 
i.e., "cut2" which search for a bump in the light curve and "cut3" in which microlensing 
fitting is performed. As a result, we find 28 microlensing event candidates from these light 
curves. 



4.1. Event Selection 

We apply a combination of simple selection criteria to reject the various spurious detec- 
tions and variable stars. These criteria are chosen empirically. 

4.1.1. cutl 

We have applied a simple selection criterion termed "point cut" for the initial variable 
objects detected in § 3.1 by using the number of detections of each object in a time series of 
frames. 

There are 3 type of events in our data set. The first one is the event in which the source 
star is not magnified in the reference frame. The second one is that in which the source 
star is magnified to its maximum amount relative to the reference image. The third one is 
in-between, i.e., the source star is magnified but not at the peak magnification relative to 
the reference image. For type 1, objects are detected only when the source is magnified as a 
positive excess. For type 2, objects are detected in almost all frames except around the peak 
as a negative excess. And for type 3, the same excess as for types 1 and 2, or the combination 
of a cluster of positive excesses at the peak and a negative base line are detected depending 
on the phase of magnification in the reference image. 

First we count the number of clusters of positive and negative detections Ncius,p and 
Ncius,n which are the consecutive detections whose separation are less than 4 observation 
frames. 
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For type 1 which has the no negative detections Ndet,n — 0, the positive detection 
should be larger than 2, i.e., Ndet,p > 2, which rejects many spurious detections and low S/N 
events. Ncius,p should be less than 4 to reject short period variables, and Ndet,p should not 
be equal to Ncius,p to reject objects with noisy sparse detections. We divide all frames into 
three regions of observation time, and calculate the detection ratio Ratioi^2,3 which is the 
number of detections Ndet,i,2,^ out of all observation frames A^/rame,i,2,3- We require that the 
Ratioi < 0.1 or Ratios < 0.1 to reject some of the long or middle period variables. 

For type 2, which has the no positive detection Ndet,p = 0, the negative detection should 
be larger than a quarter of all observation frames, i.e., Ndet.n > N frame/4: because these 
negative detections represent the base line. Ncius,n should be less than 5 to reject short 
period variables because the base hne should be stable. 

For type 3 which has both positive and negative detections, Ncius,p should be less than 4, 
and Ncius,n should be less than 5 to reject short period variables, the same as for types 1 and 2. 
Either of the criteria for typel Ndet,p > 2 and Ndet,p ^ Ncius,p, or for type2 Ndet,n > N frame/ 4 
should be satisfied. We hst all criteria of cutl in Table 1. 

These cuts reduce the variable objects from thousands (2, 000 ~ 10, 000) to hundreds 
(200 ~ 1,000), depending on each field, without rejecting the microlensing-like light curves. 
For these objects the PSF photometry is performed to construct light curves. 



The second cut "cut2" have been applied to the light curves of objects which passed 
"cutl". This cut search for a microlensing-like "bump" on the stable baseline in the light 
curve. 

First, we require that the number of photometric data points Ntotai should be > 70 
points. Next we set the time interval of 120 days as a "window" in the light curve, where 
the number of data points in the window Nin and outside the window Nout is required to be 
more than 3 and 9 points respectively. In this window, we count the number of peaks Npeak, 
which is defined to be that consecutive excess whose significance is larger than 2.5a and at 
least the significance of 2 points of these should larger than 4(7. Here the significance of each 
photometric point is calculated as follows 



4.1.2. cut2 




(13) 
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where fi and af^i are the flux and error of the ith data point, and fmed,out and a out are the 
median and variance of data points outside the window. 

Npeak is required to satisfy 1 < Np^ak < 3, and the maximum of the sum of the sig- 
nificance of the points in each peak ^ . ^^^^^ Ui should be larger than 20. We define a high 
excess as one or more consecutive data points whose significance is larger than 2.5(7, i.e., all 
excesses including "peak" defined above and one which is not significant to categorize as a 
"peak". The number of the high excesses N^i in the window should be less than 6, which 
rejects short period variables and noisy light curves. We require that the reduced chi-square 
of data outside the window x^^^ should be less than 4. But if the ratio of chi-squarc inside 
and outside the window xln/xlut larger than 15, the cut with xlut not required, so as to 
allow the high S/N long duration events to pass. We hst all these criteria of cut2 in Table 2. 

In this "cut2", 24, 543 light curves have been reduced to 1, 014 in all our GB fields. 

4.1.3. cuts 

For the light curves which passed cut2, the microlensing model fitting with equation 
(1) have been applied in cut3. In DoPHOT analysis the fittings are performed with the 
background blending flux F;,, i.e., F{t) = FoA(t) + F^, where Fq indicates a baseline source 
flux. However, in DIA, we observe only the variation of the flux from that in the reference 
image as follows 

AF{t) = FoA{t)-FoA{tref), (14) 

where t and tref are the time when the current and the reference image are taken respectively. 
Fb is canceled out in this formula. The equation (14) with four parameters (Fo, ^o, limin, ^e) 
is fltted to differential flux AF light curves in this cutS. 

The reduced chi-squares in a microlensing flt are required to be less than 3.5 to reject 
most of LP Vs and noisy light curves, although very high S /N events are sometimes not well 
fltted by the standard microlensing model because of exotic effects such as parallax. So, 
for the events whose peak ffux is larger than 450, 000 ADU, we require < 100 instead of 
< 3.5. The event ngbl-2-2717 falls into this category because of the parallax effect. After 
these cuts, 75 light curves still remain, most of which have a clear single peak and stable 
baseline. These are microlensing candidates, dwarf novae (DN) and low S/N faint LPVs. 

The main background in this event selection are DN which can be well fltted by mi- 
crolensing in the case of poor sampling. In our light curves which passed the cutl, 20 light 
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curves are clearly identified as DN. These have single or sometimes multiple asymmetric 
flares in the light curve, which usually rise quickly and fade slowly. 

To check whether an object is DN or not, we cross-referenced to the existing light curves 
in the MACHO group's database (http://wwwmacho.mcmaster.ca/) whose fields overlap 
some of ours, and found other fiares in the light curves of 5 objects, i.e., these are clearly 
DN. The DN are much bluer during the outburst. However we usually observe only with 
the red filter to increase the sampling rate, and only rarely observe with the blue (only ~ 30 
frames for each field). For 23 out of these 75 candidates, we could measure the color during 
the fiare. Five of these 23 candidates arc categorized as DN because they are very blue 
{V — I < 1.0). In 15 light curves, including 5 of the above 10 DN, asymmetric fiares are 
clearly seen. In total 20 objects are categorized as DN. 

Two of these DN are rejected at cut2. And one DN is removed at by cut3. To reject 
other DN, we are unable to use a color cut because we have no color information in most of 
the fight curves. On the other hand we have frequent sampling data points for these objects, 
which makes it easier to identify a DN from the shape. And we find the following cuts would 
reject almost all DN in addition to LPVs. 

The minimum impact parameter Mmin should be less than 1. 15 DN are fitted with a 
large Wmin- Some of LPVs are also cut here. Furthermore we impose the condition that the 
timescale should be 0.3 < is < 200 days. Two further DN are rejected here, which are fitted 
to very large event timescales (^e > 1000 days). In consequence, all DN which have been 
clearly identified are rejected by these cuts. Though 5 of these have been identified as DN 
not by a clear shape of the fiare but by their color or by the existence of other fiares in past 
data, which can be done randomly, these are rejected by our cuts. These results give us the 
confidence that all DN are rejected by these cuts. Even if a few DN are not rejected, it should 
not be a significant fraction in our results. At the same time 9 microlensing-like light curves, 
which could not be clearly identified as cither DN or microlensing events, are rejected with 
this cut. In microlensing events measured only around the peak or with low S/N, and 
i^min would bc degenerate (Han 2000). So some real microlensing events might be removed 
at this cut. These effects are also seen for the artificial events with faint source stars in our 
simulation (see § 6.3). If both the colors are taken simultaneously or the catalogue of DN is 
used, we can easily distinguish real microlensing events from DN and issue the alerts. Such 
events are usually high magnification events because their source star is very faint. If the 
baseline fiux is measured by follow up observations with larger high resolution telescope such 
as HST or VLT, the timescale is well constrained. Other LPVs are also cut here. 

31 objects have passed these criteria. In these light curves 3 low S/N LPVs still remain. 
We have rejected these LPVs directly as we are doing in real-time analysis (Bond et al. 
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2001) instead of imposing more complicated criteria. In real-time analysis we have made a 
list of variables and those objects which cross-referenced to these variables are rejected. This 
treatment does not affect the optical depth estimation because that is position-dependent 
and these positions are rejected in following analysis. We list all the criteria of cut3 in Table 
3. 

As a result, we find 28 microlcnsing candidates in our GB database during 2000. We 
confirmed that these 28 candidates didn't have any significant variability during the following 
two seasons (in 2001 and 2002). Three out of 20 candidates reported in real-time analysis 
(Bond et al. 2001) have failed in this off-line analysis because the arc not well constrained. 
And 11 new candidates are found because we have changed the threshold to detect variable 
objects on the subtracted images in §3.1. We summarize the event selection processes in 
Figure 7. 

We show light curves of these 28 candidates in Figure 8 ~ 10, where the AF data 
points have been converted to amphfications using the fitted parameters. The gap around 
JD = 2451750 in the fight curves is due to the mal-functioning of our camera system for 
~ 40 days. 

We list the positions of all candidates with ID in this analysis in Table 4. The ID 
in real-time analysis and alert ID reported in Bond et al. (2001) are also tabulated. We 
list the best fit microlensing parameters and la lower and upper limits in Table 5. The 
/-band baseline magnitude of the source star Iq is de-reddened to match the HST field by 
using the /-band extinction Aj map of each field. These fitted parameters are not biased by 
the blending effects due to nearby stars which appear in DoPHOT-type analyses (Udalski 
et al. 1994; Alcock et al. 1997a,b)and would make shorter. In DIA all other blending 
components could be subtracted. 

The light curve of the event ngbl-2-4925 exhibits an asymmetric profile due to the 
effect of parallax. The best fit parameters in fitting with parallax microlensing models for 
this event are presented in Bond et al. (2001). We will use these values for the following 
optical depth estimation. 

5. Extinction Map 

As is well known, the extinction due to dust is very significant towards the GB. This 
effect can be seen in our Color Magnitude Diagram (CMD) of the field near Baade's Window 
(left panel) and another high extinction field (middle panel) in Figure 11. The CMD is more 
scattered in the high extinction field than in the Baade's Window field due to the extinction 
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and the reddening. Information on the extinction in each region is needed to estimate the 
number of stars in this region by using the luminosity function. 

We apply Paczyhski et al. (1999) 's method which use red clump-dominated parts of 
the CMDs for determining the offsets caused by differential extinction. They have made a 
reddening map for their fields because determining the reddening Ey-i (horizontal shift in 
CMD) is easier than Ai (vertical shift in CMD). Then the extinction map arc calculated by 
the following formula (Wozniak & Stanek 1996; Ng et al. 1996; Stanek 1996): 

Ai = 1.5 X Ev-i. (15) 

Wc have made CMDs of / vs (V — I) for each of 8 x 16 sub-divided regions (3.45' x 3.45') 
by using our star catalogue. A sample of these CMDs of different subregions in the image near 
Baade's window is shown in the right panel of Figure 11. In this figure, differential reddening 
is clearly seen. We select the red-clump-giant region in the CMD following Paczyhski et al. 
(1999) as indicated in Figure 11. The parallel lines whose slopes are the reddening vectors 
(Stanek 1996) are given by 

7- 1.5(y- 7) = 1.2,13.25. (16) 

We estimate the mean vahics of (V — 7) and 7 for red clump giants in each subregion, 
which correspond to the correlation between differential reddening E{V — I) and differential 
extinction Aj. These values in a field are plotted in Figure 12, where the best fitted line is 
given by 

7= (1.45±0.12)(l^-7) + 12.7. (17) 

This slope is consistent with Stanek (1996). The mean colors {V — I)mean is more reliable 
than 7. So we estimate Ey-i from {V — I) as follows. Wc compare (V — I)mean from our 
data and the reddening Ey-i map (centered on (18:03:20.9, —30:02:06)) calculated by Stanek 
(1996) in the overlap regions (ngb2-l, ngb2-2 and ngb3-3), where the zero point is based on 
the determination by Gould et al. (1998) and Alcock et al. (1998). A correlation between 
{V — I)mean Ey-i is shown in Figure 13 and equation (18). 



Ey_j = {V-I)mean-i-i8± 0.05. (18) 

Using equations (15) and (18), we have made an 7-band extinction Aj map from the 



-16- 



{V — I)mean map in our data as shown in Figure 14. The error in Ai is about 0.07 mag. This 
map is used to make a luminosity function for each region. 

6. Optical Depth 

To determine the microlensing optical depth from detected events we estimate our de- 
tection efficiency by Monte Carlo simulations. Here we add artificial microlensed star images 
to the real subtracted images instead of a full Monte Carlo approach. Then we count the 
number of events which passed the same event selection applied in the real event selection. 
Then we estimate the optical depth towards the Galactic bulge. 

6.1. Luminosity Fiinctions 

To add the artificial microlensing events, we need the Luminosity Functions (LF) of 
source stars in GB fields. We use the deepest observations with the Hubble Space Telescope 
(HST) from Holtzman (1998), which measured stars in Baade's window down to / (F814W) 
~ 24. 

A completeness correction has been applied for the HST LF for fainter stars. The 
bright end of the HST LF is poorly defined because of the small field of view (~ 5 arcmin^). 
Meanwhile the MOA data have a lot of bright stars because of the wide field of view (~ 1.2 
deg^), but are poor for fainter stars. So we combine the HST LF and MOA LF at / = 15 ~ 16 
as shown in Figure 15. Here the MOA LF is made from stars in field ngb2-2 (0.4 deg^) which 
includes the HST field. The offsets in the photometry of stars in each subregion due to 
differential extinction are corrected to match that of the HST field {Aj = 0.742 mag). We 
use the MOA LF for 7 < 16 and the HST LF for 7 > 16 in this analysis. We use this 
composite LF for all our GB fields assuming that the morphology of the LF varies little 
between the small HST field and the large MOA fields. 

6.2. Simulation 

In each of our 42 (14 fields x 3 chips) subfields, the density of stars, samphng rate, and 
observational conditions are different. So we estimate the detection efficiencies individually 
for each subfield by a Monte Carlo simulation. We generate 429,000 artificial microlensing 
events in each subfield for this purpose. 
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A given event has 4 parameters: (limin, ^e, ^o, -^o)- We generated artificial events with the 
timescale = (0.3, 0.5, 1, 3, 5, 10, 20, 40, 60, 80, 100, 150, 200) days. For each tg, wc generated 
events with source extinction-corrected /-band magnitude between 12.0 < Iq < 22.9 mag at 
intervals of 0.1 mag uniformly (the corresponding flux values Fq were estimated by taking 
the extinction at each position into account). For each value of and /q, 300 events were 
generated (i.e. 429,000 events for each subfield). 

Those parameters which follow the well known distribution such as the position (x, y) , 
the peak time to and the impact parameter -Umin are selected at random uniformly in the CCD 
chip (36 < .X < 2042 pixels and 5 < a; < 4090 pixels), in the observation period (2451645.1 < 
JD < 2451863.9) and in < u^in < Uth, respectively. Here, Uth = 10-°-'^^^°-^'''^ {Iq > hh) and 
Wth = 1 (/q < Ah). We set Ith — 16 mag (for ngbl,4,5 and 10), 16.5 mag (for ngb2,3,6,8 and 
11) and 17 mag (for ngb7,9 and 12) depending on the mean extinction in their fields. We 
ensured that the events with Umin > Uth are not detectable in our experiments from former 
simulations. This Uth is introduced to enlarge the statistics for the faint source events. The 
slope of 1/iith is roughly consistent with the slope of our combined GB LF, ip{I), which 
means the number of generated events is roughly following the LF for faint source events. 

Many other components which are expected to affect the event detection, e.g., bad pixels, 

seeing, extinction, sample rate, star density, focusing, tracking, sky background, clouds, PSF 
variation in the focal plane, differential refraction, variable stars, saturated stars, cosmic 
rays, satellite, asteroids, systematic residuals and any other unknown systematic noise, are 
automatically simulated because we use real subtracted images in this simulation. 

To generate artificial microlcnsing events, we add the artificial differential star images 
on a time series of real subtracted images. The differential stars are made from the local 
PSF which was formed by the same procedure outlined in § 3.2, and scaled by the differential 
flux expressed in equation (14). 

In generating the artiflcial image, the photon noise is taken into account following Pois- 
son statistics with an additional flux in the current image. The ADU flux of the artiflcial 
star in the current image FoA(t) in equation (14) is calculated in the scale of the reference 
image. So, to estimate the photon noise from the additional flux, this additional flux should 
be transformed to the number of electrons in the current image given by FoOq/G, where G 
and oo are the gain in ADU/e~ and a scale factor in the kernel respectively. Then the elec- 
tron signals are simulated randomly following a Poisson distribution with the mean value of 
FoOo/G. The resultant simulated random electron signals are rescaled back to the reference 
image in ADU by multiplying by G/oq. 

The photon noise from the sky and the blending stars in the current image and reference 
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image are already included in the real subtracted images. Furthermore the systematic noise 
from the poor subtraction due to the low S/N is also included. These noise sources heavily 
depend on the sky condition, seeing, tracking of the telescope and star density of each frame. 
So it is complicated to simulate with full artificial images. 

We cut out a small sub-image of 23 x 23 pixels from the subtracted image at a randomly 
selected position, and put the corresponding differential flux AF{t) onto each sub-image. A 
sample of the cut raw subtracted image (left panel) and the generated artificial differential 
image (right panel) are shown in Figure 16. 

For a series of generated artificial images, we apply the same variable object detection 
process as in §3.1, and the same PSF photometry as in §3.2. Saturated pixels (> 30,000 
ADU) are rejected as in the real analysis. The resultant light curves are stored in the 
database of artificial events. Sample artificial light curves are shown in Figure 17. 

We count the number of artificial events n{tE, lo) which passed all criteria for the real 
event selection for each and Iq. Then we calculate our detection efficiencies. 



6.3. Detection EfHciencies 

We show our detection efficiencies as a function of the event timescale, the source mag- 
nitude and the minimum impact parameter. 

6.3.1. Detection Efficiency as a function of Event Timescale 

The detection efficiency £o(^e, -^o) for each sample with and Jq is given by £o(^e, -^o) = 
n{tE, /o)/300, where n{tE, Iq) is the number of artificial events which passed our event selec- 
tions. Then efficiency e{tE) for tE is integrated by weighting by the LF ■0(7) as follows, 

<tE) = ^ eo{tE,lMl)ut^dI, (19) 

^^23 J 1=12 

where 

/•/=23 

A^23= / i'{l)dl. (20) 

Jl=12 

We show one of the resultant detection efficiencies as a function of for bright stars 
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(7 ~ 14, dashed line) and for all possible sources {I < 23, solid line) in Figure 18. Most 
events with bright source stars whose timescale is longer than 3 days can usually be detected. 

On the other hand, the total event detection efficiencies (sohd line, I < 23) are very low. 
This is because most source stars are fainter than our observational limiting magnitude of 
I ~ 17. These faint source stars can be detected by microlensing only when they are highly 
magnified. The possibility that such a high magnification (i.e. small impact parameter limin) 
event occurs is small. Furthermore, such faint source events cannot be detected any more 
for shorter timescale events. 

6.3.2. Detection Efficiency as a function of source magnitude 

We present the detection efficiencies averaged over all subfields as a function of the source 
/-band magnitude Iq for various event timescales in Figure 19, where /-band magnitudes are 
extinction corrected to match those of the HST field {Aj = 0.742). The efficiencies are very 
high for brighter source events, as expected. One can see the fall off at the bright end of 
this efficiency curve due to the saturation of the source star images. The fact that this effect 
appears in our simulation is one of the pieces of evidence that this simulation is realistic. 
On the other hand efficiencies fall off for fainter stars as expected. The maximum of our 
detection efficiency is not high, even for bright source events, because our data have a gap 
for ~ 40 days due to the mal-functioning of our camera system. 

Comparing each of the timescales, for the long timescale events, the efficiency falls 
off more significantly than that of the shorter timescale events at the bright end of these 
curves. This is because the number of saturated points is large and we required a stable 
flat baseline to identify the event. In the long timescale events (^e > 50 days), the stable 
flat baseline in the light curve is very short. These effects are significant for brighter source 
and longer timescale events and vice versa. Meanwhile, at the faint end, the efficiencies are 
inverted compared with those at the bright end. For long timescale events the duration with 
a significant excess in the light curve is longer than that of shorter timescale events which 
cannot be detected any more. We note there is a small but still non- negligible efficiency 
down to / ~ 22. This is the main improvement in the analysis by changing from DoPHOT 
to DIA. 

To emphasize this effect, we show the relative expected event rate as a function of Iq for 
various event timescales in the top-left panel of Figure 20. These are estimated by multiplying 
Figure 19 by the combined GB LF '0(/) (Figure 15), i.e., the integrand of equation (19). The 
expected source distributions are peaked around / ~ 19 and go down to / ~ 22, although the 
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observational limiting magnitude of our telescope and camera is / ~ 17. This does show the 
dramatic capability of DIA. Especially note that this effect is increased by the high sampling 
rate strategy that we have adopted. Though the sensitivity of the MACHO group (Alcock et 
al. 2000) became 2 magnitudes deeper than their observation limiting magnitude of V ~ 21 
by DIA, that of MOA is ~ 5 magnitude deeper than our observation limiting magnitude of 
/ ~ 17. (Of course, this effect also depends on the shape of the luminosity function around 
the observation limiting magnitude for each experiment.) The longer the timescale is, the 
fainter the source events that could be detected. 

In the top-right panel in Figure 20, the observed /-band baseline magnitude Iq distribu- 
tion is shown as the histogram with corresponding expected event rates scaled to match the 
histogram. Here the observed Jq is de-reddened to match that in the HST field. The expected 
and observed distributions are in good agreement for the observed range (5 < < 100 
days) . 

As mentioned in § 4.1.3, many possible microlensing events are rejected by the condition 
< 200 days in cut3. These events have a significant increase in A F, but their parameters 
could not be well constrained, because the light curve have been measured only around the 
peak. If the baseline flux is measured by the follow-up observations, t-E could be constrained. 
If we make a catalogue of all variable stars and DN as now we arc doing, or if we observe 
with two colors, wc could easily identify the real microlensing events and issue alerts without 
the < 200 days criterion. 

We show the expected event rates without the cut t^, < 200 days in the bottom-left 
panel of Figure 20. In this figure we can see a significant increase in the expected event rate 
for the longer timescale events with the dimmer source stars which intrinsically should have 
high amplification. We can see this effect in the minimum impact parameter distribution 
in the following §6.3.3. On the other hand, no significant change can be seen for the short 
timescale events, though the duration of the excess in the light curve of these events are 
as short as that of the long timescale faint source event. So, the efficiency for the short 
timescale events is not affected by this criterion. 

To compare these results with the traditional DoPHOT analysis in which the source 
star should be resolved in the reference image, we show the same distributions by using 
only the MOA luminosity function (filled circle points in Figure 15) as the source luminosity 
distribution in the bottom-right panel of Figure 20. This is not actually the same case as with 
DoPHOT, because the DIA can improve the photometric accuracy, even for stars resolved 
in the reference image. Furthermore there are many faint source events which are blended 
by bright resolved stars. These effects will cancel each other. In any case they are relatively 
small effects when compared with the differences in the DIA and DoPHOT analysis. So 
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this distribution should not be much different from the actual one and make a meaningful 
comparison. In this figure, the difference from DoPHOT (bottom-right) to DIA (top-left) is 
clear. The peak of this distribution is at / ~ 15, which corresponds to the second peak in 
the top-left panel of Figure 20. The DIA method increases the number of events by more 
than a factor of 2. In particular we can detect many events in which the source stars are 
fainter than the observational limiting magnitude by using DIA. 

6.3.3. Minimum Impact Parameter Distribution. 

Because of the simple geometry of microlensing, the probability distribution of the 
minimum impact parameter Umin is uniform. In actual experiments, however, the distribution 
of the detection efficiency in Mmin, i-c. the expected event rate distribution, is biased towards 
smaller 'Umin, because the high amplification makes it easier to detect the event. 

This effect can be seen in our observed u^ain distribution (histogram) and the estimated 
one in the simulation (lines) for various timescales which are scaled to match the histogram 
at Minin = 0.05 as shown in Figure 21. Here the lines are the mean efficiencies for all our 
fields. These observed and estimated distributions are in good agreement for the range of 
the observed timescales (between 5 ~ 100 days). 

In this figure we can see that half of the detected events have a small minimum impact 
parameter (wmin < 0.1), though the fraction of events with u < 0.1 is only 10% in the 
DoPHOT analysis (Alcock et al. 1997a) and 30% in the recent DIA analysis (Alcock et al. 
2000) by the MACHO group. This is because our sampling rate is higher (5 ~ 6 times day~^) 
than theirs (once per day). 

If events which fail by the criterion t^ < 200 days, are included as mentioned in the 
previous section, then 5 ~ 15% more high magnification events would be expected. We show 
the estimated limin distribution for the typical timescale of = 40 days without (solid line) 
and with (dashed fine) this < 200 days cut in Figure 22. And we also show the same 
distribution with tE < 200 days cut for brighter source stars (/ < 17), which correspond to 
the DoPHOT analysis. 

From this figure, it is clear that our analysis with DIA (solid line) can detect the 
high magnification events more efficiently than a DoPHOT- type analysis (dot-dashed line). 
Furthermore it will be possible to detect more (~ 10%) high magnification events without 
the tE < 200 days cut (dashed fine) in the near future. 
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6.4. Optical Depth 

Here we estimate the optical depth towards the Galactic bulge by using the observed 
events and our detection efficiencies. 



6.4-1- Optical Depth Estimation 

The optical depth r is defined as the probability that any given star is microlensed with 
impact parameter limin < 1 at any given time. The r can be estimated by 

_ tE,i 

^ ~ 2N,T, ^ e{tE,i) 

where Ng is the total number of source stars in the observation fields and Tq is the duration 
of observation of this analysis in days. tE,i is the event timescale for the ith event and £:(tE,i) 
is the detection efficiency for a given timescale. 

We estimate Ns for / < 23 by using the number of bright stars (/ < 15 which is 
extinction corrected) whose observation efficiency is nearly 100%, and our combined LF. 
This is because our detection efficiencies have been estimated by taking all stars with 7 < 23 
into account. 

We estimate the optical depth by using the observed timescales with Ng ~ 250 million 
stars (/ < 23), Tq = 219 days for ngbl ~ 12. We do not include fields ngbl3 and 14 in this 
analysis because these two fields are far from Baade's Window. For events with timescales 
within 0.3 < tE < 200 days, we estimate the optical depth as 

T = 2.63i°:5l X 10-^ (22) 

where the lower and upper limit of this value are estimated by a Monte Carlo simulation 
following Alcock et al. (1997a,b, 2000) and assuming a Poisson distribution. We have simu- 
lated the 100,000 "experiments" for each A^exp which is the number of expected events and 
chosen at 0.5 intervals between and 60 events. In each experiment the number of detected 
events N are selected by following the Poisson statistics with the mean value of Nexp- For 
each simulated event, we randomly selected one of our observed event timescales. We have 
estimated the probability P{t{N) > Tobs) that the optical depth t{N) in each simulated 
experiment is larger than the observed one Tob^, and the mean optical depth {T{Nexp)) for 
each distribution with Ng^p. We show this probability distribution in Figure 23, from which 
we have estimated the la confidence limit of the optical depth. 



(21) 



-23- 



This estimate gives errors which are a lower hmit, because errors in tE are ignored. 
We found a large uncertainty and systematic bias in the "INPUT" (hereafter tEin) and 
measured "OUTPUT" (hereafter tEout) in our artificial events generated in § 6.2. We show 
the relation between tEin and the mean value of tEout for various source magnitude in Fig. 
24. As we can see in Fig. 24, the measured tEout tend to be larger in small tEm and smaller 
in large tEin- This is because we cut events with tEin outside 0.3 < tE < 200 days. This effect 
is larger especially for fainter source events because of their large scatter in tEout- However, 
for relatively brighter source events with 10 < tsin < 60 days, there is only a quite small 
bias that tsout tend to be shghtly larger. 

We carried out a full Monte- Carlo simulation to measure the error in the optical depth 
and bias by using these artificial events. We used lists of measured OUTPUT tEout of 
artificial events for each "INPUT" value of tEin and /q, which consist of about a hundred 
thousand of artificial events for each subfield. The number of expected events for each Iq 
is proportional to ~ £o(^e, -^o)V'(-^o)^*th, i-e-, the integrand of equation (19). Following this 
number of expected events for each Jq, we picked up OUTPUT tEout randomly from the lists 
and put them into the expected OUTPUT tEout distribution -Dout(^Ein) until the total number 
of events become ~ 40, 000 for each INPUT tEin- We show a example of this distribution for 
tEin = 40 days, i.e. Dout(40), in Fig. 25. 

We used our observed tE distribution as the INPUT tEin distribution for the simulation. 
We don't have any artificial events with exactly the same tE as that of observed events. So 
we linked each event to artificial events which have a similar INPUT tEin value. For example, 
event ngbl-2-2745 (tE = 41.3 day) was linked to artificial events with tEin = 40 day, and 
ngbl-2-5076 (tE = 51.8 day) was linked to artificial events with tEin = 40 and 60 days with 
50% probabihty. The linked INPUT tEin of artificial events were listed in the INPUT tEin 
distribution Di^. 

We performed the same simulation as above by using Din and -Dout(^Ein), i-e-, we sim- 
ulated 100,000 "experiments" for each A'exp- For each simulated event, we randomly select 
one of tEin from Din and choose one of the tsout from corresponding OUTPUT distribution 
Dout(iEin) at random. Then we measure the INPUT and OUTPUT optical depths (Tin and 
Tout) from tEin and tEout respectively, for each experiment. We ensured that the Tin in the 
case of iVexp = 28 is the same as the observed value Tobs = 2.63 x 10~^. This is expected 
because Di^ was made from our observed tE distribution. 

We plot the difference between the mean value of Tin and Tout in Fig- 26. We can see the 

existence of the small bias between them. The Tout tend to be slightly larger than Tin. From 
this relation, the real optical depth is estimated to be t = 2.59 x 10~^ from the observed 
optical depth of Tobs = 2.63 x 10~^. 
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To estimate the error in this optical depth, we have used the same method as above 
but with Tout instead of Tin- The standard deviation of Tout is ~ 16% larger than that of 
Tin for each iVexp- We have estimated the probability P(rout(-^) > Tobs) that the optical 
depth rout(-^) in each simulated experiment is larger than the observed one Tobsj and the 
mean INPUT optical depth < rin(iVexp) > for each distribution with N^xp- We show this 
probability distribution in Fig. 27. From this we have estimated the la confidence limit of 
the "INPUT" optical depth Tin. These la errors are about 15% larger than errors from the 
method with Poisson statistics only. The resultant optical depth and errors are 



T = 2.59l°:^t X 10-^ (23) 

This T is mostly due to disk and bulge (bar) stars. The optical depth in the direction of 
the Galactic bulge due to halo objects of any kind is only ~ 0.13 x 10~^ (Griest et al. 1991). 
This T is underestimated because some fraction of the source stars are foreground disk stars, 
for which the optical depth is considerably lower. 



6.4-2. Disk Contribution 

The optical depth estimated in the previous section is underestimated. In the number 
of stars in our fields A^^j some fraction of the stars are probably foreground Galactic disk 
stars. The optical depth for the disk source stars (so-called disk-disk events) is quite small. 
There are also background disk stars which would have a higher optical depth, but our line 
of sight towards the bulge is several hundred parsecs out of the Galactic plane on the far side 
of the bulge. So most of the disk contamination is from foreground disk stars. The fraction 
of disk stars in our fields is rather uncertain. 

We estimate the fraction of disk stars out of all stars, /disk in our fields. We use the 
non-rotating triaxial bar models with the bar inclination angle of ^ = 20° whose density 
profile as a function of the distance D from the Sun is given by (Han & Gould 1995a; Alcock 
et al. 2000) 



M / U)2 

^^^"^ = 20:65^^"n"T 
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Here the coordinates {x', y') are measured along the longest axis and another axis of the 
bar in the Galactic plane. The x'-axis is aligned at an angle 9 from the line of sight to the 
Galactic center from the Sun, with the near side of the bar in the positive-/ quadrant. The 
^-axis is as usual the height above the Galactic plane. These Galactocentric coordinates 
(x', I/', z) are given by x' = Rq cos 6 — D cos h cos(/ + 6'), y' = Rq sin 9 — D cos b sin(/ + 9) and 
z — D sin b where Rq is the distance to the Galactic center from the Sun taken to 8.5 kpc. 
a = 1580 pc, b = 620 pc and c = 430 pc define the bar scale lengths, and M = 1.8 x 10^° Mq 
is the total bar mass. 

Wc use the standard double-exponential disk whose density profile is given by (Han & 
Gould 1995a; Alcock et al. 2000) 



where R — Rl + D"^ cos^ h — 2RqD cos h cos I and z — Dsmh are the disk cyhndrical Galac- 
tocentric coordinates, D is the distance from the Sun and {I, h) are the Galactic latitude 
and longitude, hz — 325 kpc and Rd — 3.5 kpc are the disk scale height and length, 
pdo = O.O6M0 pc~^ is the density constant chosen to match the density in the solar neigh- 
borhood. 

These models give optical depths towards Baade's window {I, h) — (1°, —3°. 9) of 1.2 x 
10~^ from the bar, and 0.6 x 10~^ from the disk. By using these density models and Kiraga 

& Paczynski (1994)'s luminosity function with /9 = — 1 (see § 7), /disk is estimated as ~ 23%. 
This value is consistent with the value used by the MACHO group of fdisk ~ 20% (Alcock et 
al. 1997a) and fdisk ~ 25% (Alcock et al. 2000), where most of our GB fields are overlapping 
with their fields. This gives 



This optical depth is consistent with the previous observations of 3.3li;2 x 10 ^ from 9 events 
by DoPHOT analysis (Udalski et al. 1994), 3.9]'^}'2 x 10~^ from 13 events in the clump giant 
subsample by DoPHOT (Alcock et al. 1997a) and ?>.2?,tQ^^ x 10"^ from 99 events by DIA 
(Alcock et al. 2000), while slightly higher than 2.0121 x 10"^ (Popowski et al. 2000) and 
2.23l2J^ X 10"^ (Popowski 2002). 

This measured optical depth must be regarded as a lower limit of the true value because 
our observations are only sensitive to the events with 0.3 <t-^ < 200 days. Nevertheless this 
value is still higher than those predicted by most Galactic models whose mass and inclination 
of the bar are consistent with other observations. We note that our observed optical depth 




(26) 



- 3 36+^-^^ 



xl0-6[0.77/(l-/disk)]. 



(27) 



-26- 



values are averaged over 16 deg^ (12 fields) around Baade's Window (Z, 6) = (1°.0, — 3°.9), 
while in most calculations the optical depth is estimated exactly towards Baade's Window. 
However we ignore this difference in this paper since it is negligible in most models. Even 
the smallest inclination angle and a large bar mass have been reported to be insufficient to 
produce an optical depth greater than ~ 2.5 x 10~® (Peale 1998). 

7. Timescale Distribution 

The observed timescale distribution depends on the mass function and the velocity 
dispersion of the lens. To make meaningful models of Galactic structure, both the optical 
depth and the timescale distribution should be consistent with the observations. Peale 
(1998) and Mcra ct al. (1998) show that it is difficult to reproduce the observed timescale 
distribution of Alcock et al. (1997a) with the existing Galactic models. We show our observed 
timescale distribution in Figure 28. 

To be compared, we also plotted the expected timescale distribution for a fixed bar 
(equation (24)) and disk (equation (26)) density model with various mass functions in Figure 
28. These distributions arc corrected by our detection efficiencies and normalized to the 
number of observed events. Where we followed the method of Kiraga & Paczyhski (1994), 
Han & Gould (1995a) andAlcock et al. (2000) assuming the mean velocity v and variance a 
of the components in the bar and disk for each direction {y and z) as {vbar,y, Vbar,z) — (0, 0) 
and {abar,y,(^bar,z) = (110,110) for the bar lens and source, {vdisk,y,Vdisk,z) = (220,0) and 
i(7disk,y,crdisk,z) = (30,30) for the disk lens, {vo,y,Vo^^) = (220,0), {ao,y,ao,z) = (0,0) for 
the observer in kms^^ (Han & Gould 1995a; Alcock et al. 2000). For these bar and disk 
models, we evaluate timescale distributions towards Baade's window with the following five 
mass functions: (i) Scalo (1986) 's Present Day Mass Function (PDMF), a (5-function (ii) at 
M = O.IMq and (iii) at M = l.OM©, and the power-law 0(M) oc M"" (iv) with a = 2.3 
and the low mass end Mi — O.IM0, (v) with a — 2.0 and Mi — O.OIMq, where the function 
(v) represents the brown dwarf rich mass function (Alcock et al. 1997a). 

The timescale of detected events is distributed in the range 5 < < 100 days and 
centered around ~ 30 days. This feature is consistent with that from the MACHO group 
(Alcock et al. 2000) though their distribution is slightly sharper around the mean than ours. 
The distribution of Alcock et al. (2000) is well fitted by the timescale distribution expected 
from the Scalo (1986) 's PDMF except for some fraction of long timescale events, though 
those models do not explain the large observed optical depth. For our distribution in Figure 
28, Scalo (1986)'s PDMF also seems to be more reasonable than the others. Although 
our detection efficiency is sufficiently high for the short timescale events (^e ~ 0.3 days) 
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because of our frequent sampling (5 ~ 6 times day~^), there are only two short timescale 
events (^e < 4 days). The timescale distribution traces the mass distribution, for the 
most common events is ~ 7(M/O.1M0)^/^ days for both bulge and disk lenses (Han & Gould 
1995a). The number of observed short events is much smaller than that expected from the 
brown dwarf rich mass function (v), which is consistent with Alcock et al. (2000). However, 
this topic is quite complicated. The measurement of has large uncertainty in itself, and 
its distribution depends on the unknown kinematics ofthe sources and lenses. A detailed 
analysis of this topic is beyond the scope of the present study. The identification of brown 
dwarfs may be possible only with larger statistics of higher S/N events in the future. 

We can also see some fraction of long timescale events (^e > 70 days) as reported in 
Alcock et al. (2000), which cannot be reproduced with Scalo (1986) 's mass function (i) and 
the contemporary Galactic models (Han & Gould 1996). Such long timescale events could 
be produced when the lenses are heavy, moving at low transverse velocity, or in the middle 
of the line from observer to source. If both the lens and source are in the disk, so called 
disk-disk lensing, then the timescale would be long. The probability of this is constrained by 
star counts, and it is very small. There might be some unknown population of dynamically 
cold or massive dark objects, such as white dwarfs or neutron stars in the Galactic disk or 
bulge. More discussion on this can be seen in Alcock et al. (2000). 

We show our observed timescale distribution, corrected by the detection efficiency, i.e., 
the expected true t-^ distribution in Figure 29. This distribution is similar to that in Alcock et 
al. (2000) except the part for the short timescale events. Though their distribution is sharply 
truncated at ^e ~ 4 days except for one short event with ^e = 1-4 days, our distribution is 
flat to 2 days. However this difference is not significant because these shorter two bins are 
based on only two events, and the amount is still small. 

The contribution to the total optical depth of the observed timescale distribution is 
given in Figure 30. The contributions of the short timescale events to the total optical depth 
are quite small. 

In any case, we need more observations to investigate the mass function and Galactic 
structure in more detail. 



8. Discussion and Conclusion 

We have re-analyzed the sample of subtracted images that were derived from the real- 
time DIA of GB observations obtained by MOA during 2000 (Bond et al. (2001)). In this 
analysis we have found 28 microlensing event candidates in our 12 GB fields. The DIA is 
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more suitable than DoPHOT analysis for our purpose, since the former method can detect 
the luminosity variation at any position, even where no star was previously identified. 

We have used a Monte Carlo simulation to estimate our event detection efficiencies. 
By using these efficiencies and timescales of our 28 detected events, we have estimated the 
optical depth towards the GB for events with timescales within the range 0.3 < < 200 
days as Tq — 2.63^g'gg X 10 ^, where the statistical uncertainty of the optical depth have 
been estimated by Monte Carlo simulation only using the Poisson statistics. This r is 
overestimated and these errors in r are underestimated. We found the systematic bias 
and large uncertainty in the OUTPUT tEout in our simulation. By taking these bias and 
uncertainty in tsout into account, we get 



= 2.59t'oit X 10-^ 

The GB microlensing optical depth, in which the disk source stars component (/disk = 23%) 
is taken into account, is given by 



Tbuige = 3.36l^ii X 10-6[0.77/(l - /disk)]. 

This value is consistent with the previous observations of 3.3^i;2 x 10~^ from 9 events by 
DoPHOT analysis (Udalski et al. 1994), 3.9;?:2 x lO'^ from 13 events in a clump giant 
subsample by DoPHOT (Alcock et al. 1997a) and 3.23;o:52 x 10"^ from 99 events by DIA 
(Alcock et al. 2000), while slightly higher than 2.01:^:4 x 10"^ (Popowski et al. 2000) and 
2.231^:1 X 10"^ (Popowski 2002). 

This observed r must be regarded as a lower limit of the true value, since our observations 
are only sensitive to events with 0.3 < < 200 days. Nevertheless this value is still higher 
than those predicted by most Galactic models whose mass and inclination of the bar are 
consistent with other observations. Even the smallest inclination angle and a large bar mass 
could not reproduce r greater than ~ 2.5 x 10"^ (Peale 1998). In Evans & Belokurov (2002) 
Freudenreich's model can reproduce the high optical depths toward the GB ~ 2.5 x 10~^ 
which is about 1 a level of our estimate. 

The smallest inclination produces the largest optical depth. However the bar inclination 
has been reported to be in a wide range between 10° — 45° summarized in Table 7 in Alcock 
et al. (2000). The optical depth also depends on the mass of the bulge or bar. However 
various observations provide conflicting values Miyuige\bar — 0.7 — 2.8 x 10^*^ Mq (Zhao & Mao 
(1996); Dwek et al. (1995); Holtzman (1998); Han & Gould (1995b); Blum (1995); Zhao 
& Mao (1996)). Very small inchnations {9 = 11° with Mi,ar = 2.0 x lO^^M© (Zhao & Mao 
1996) and 9 = 12° with Mi,ar = 2.5 x IO^^Mq (Gyuk & Crotts (1999))) or a very heavy 
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mass (~ 3.6 x lO^^M© with 9 = 20° (Gyuk & Crotts (1999))) are required to account for the 
observed optical depth in this analysis and in Alcock et al. (2000). Binney & Evans (2001) 
estimate the minimum total mass in baryonic matter within the Solar circle to be greater 
than ~ 3.9 x lO^^M© from r = 2.0 x 10^®, and such high baryonic contribution is consistent 
with implications from hydrodynamical modeling and the pattern speed of the Galactic bar. 

Such a high mass would imply low halo MACHO fractions (Gates ct al. 1996). A massive 
bulge puts tight constraints on the contribution of the disk to the rotation curve at small 
radii. A small disk, however, leaves more room for the halo. Since microlensing results 
towards the LMC fix the MACHO content in the halo, a massive halo implies a smaller 
MACHO fraction. 

The uncertainties in the Galactic bar orientation, the bar mass and the stellar mass 
function are still large. So the optical depth we have derived might yet be explained by 
other models. Further discussion can be seen in Alcock et al. (2000). 

We showed our observed timescale distribution which is not biased by blending. Though 
our statistics are smaller, this seems to be consistent with the one previously presented by 
the MACHO group (Alcock et al. 2000). The number of short timescale events was quite 
small, in spite of our high detection efficiency for events down to ~ 0.3 days. 

A significant number of long timescale (tg > 70 days) events have been detected, as 
reported previously by the MACHO group (Alcock et al. 2000). These could not be explained 
by any current Galactic model (Han & Gould 1996). These might be a heavier remnant 
component, such as white dwarfs, or some dynamically cold component. Either way, we 
need more observations to investigate the mass function and Galactic structure in greater 
detail. 

We have shown how efficiently MOA can detect the high magnification events in which 
the probability of detecting extrasolar planets is high and find 50 ~ 60% of all detected 
events have high magnification (limin < 0.1). This fraction is much higher than 10% from 
the DoPHOT analysis (Alcock et al. 1997a) and 30% from recent DIA analysis (Alcock et al. 
2000) by the MACHO group. This is because our sampling rate is higher (5 ~ 6 times day ~^) 
than theirs. These results support our belief that high frequency observations and analysis 
using DIA, that MOA is currently carrying out, can detect high magnification microlensing 
events very efficiently, even with a small telescope. 
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Fig. 1. — Half width of the region containing 68% of residuals which are a ratio of the actual 
scatter to the photon noise estimate as a function of SdeVpsf (filled circle) and the same plot 
with renormalization (cross) . The dashed line indicates the best fit for the residual plot with 
3-degree polynomial. 
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Fig. 2. — Distribution of residuals (upper) and renormalized residuals (lower) for 90,000 
individual measurements of constant stars as histograms. The dotted lines indicate the best 
fit Gaussian distributions centered on 0. The half widths of the regions containing 68% of 
the residuals a are 1.369 and 1.005 respectively. 



-35- 



I 1 1 1 r 



-"I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



_i I I I I I I I I I I I I I I I I I I I I I I I I I I I i_ 



13 



14 



15 



16 



17 



18 



19 



I (mag) 



Fig. 3. — Half width of the region containing 68% of the residuals which are a ratio of 
the actual scatter to the error from the photon noise (filled circle) as a function of /-band 
magnitude and the same plot with renormahzation (cross). The data in this plot come from 
90, 000 individual measurements of 1, 000 constant stars. 
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Fig. 4. — Calibration of MOA Red measurements to the standard / band magnitude for the 
three CCDs (chipl, 2 and 3 from top to bottom) with the color term. 
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Fig. 5. — Calibration of MOA Blue measurements to the standard V band magnitude scale 
of the three CCDs of the camera M0A-cam2 with the color term. 
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Fig. 6. — Calibration of MOA Red measurements to the standard / band magnitude scale 
of the three CCDs of the camera M0A-cam2 without a color term. 
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Fig. 7. — A flowchart of the procedure for selection of microlensing candidates. 
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8. — Light curves of microlensing event candidates. 
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Fig. 9.— 



Light curves of microlensing event candidates. 
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Fig. 10. — Light curves of microlensing event candidates. 
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Fig. 11. — CMD of the Baade's window field (left panel), a higher extinction field (middle 
panel) and two subregions in the Baade's window field (right panel). The extinction Aj in 
each subregion (in right panel) are 0.761 (filled circle), 1.645 (open circle). The box in the 
right panel encloses the red clump giant region defined by / — 1.5{V — I) — 1-2, 13.25. 
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Fig. 12. — Correlation between the mean of {V — I) and / for red clump giants in each 
subregion. The fitted line is / = 1.45(1^ ~ + 12.7. 
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Fig. 13. — Correlation between {V — I)mean from MOA and Ey-i from Stanek (1996) in the 
overlap regions. The fitted hne is Ey-i = (V^ — I)mean — 1-18 i 0.05. 
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Fig. 14. — /-band extinction Aj map of 14 Galactic Bulge fields. The darker regions indicate 
higher extinction, lighter regions mean lower extinction. We can see the galactic disk, in which 
the extinction is high, in the diagonal line from the top-right to bottom-left. 
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Fig. 15. — Combined luminosity function of Baade's window from MOA (filled circle) and 
HST (filled triangle). For MOA data the offsets due to differential extinction are corrected 
to match that of the HST field. We use MOA data for 7 < 16 and HST data for 7 > 16 in 
this analysis. 
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Fig. 17. — Sample light curves of artificial microlensing events with = 20 days and /q ~ 17. 
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Fig. 18. — MOA detection efficiency as a function of for tlie ngbl-2 subfield for all source 
stars (/ < 23, solid line) and for bright sources (/ = 14, dashed line). The efficiencies slightly 
differ for each subfield because of differences in the sampling rate, the star density and the 
extinction. 
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Fig. 19. — MO A detection efficiencies averaged over all subfields as a function of source 
/-band magnitude Iq for the various event timescales as indicated in Figure. 
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Fig. 20. — Expected relative event rates as a function of the source /-band magnitude Jq for 
the various event timescales as indicated in Figure. Top-left: Expected event rates in this 
analysis. Top-right: Histogram of the observed baseline magnitude Jq with corresponding 
expected event rates scaled to match the histogram. The expected and observed distributions 
are in good agreement for the observed range of (5 < < 100 days). Bottom-left: 
Expected event rates without the cut of tE < 200 days in cut3. A significant increase can be 
seen for the longer timescale events with the dimmer source stars. Bottom-right: Expected 
event rates for the case that the sources are only resolved stars in the reference image, which 
corresponds to the event rate by DoPHOT-type analysis. This is shown in comparison with 
that by DIA. 
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Fig. 21. — Observed u^i^ distribution (histogram) and the estimated distributions (hnes) 
scaled to match the histogram at Uj^in — 0.05 for various timescales as indicated in the figure. 
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Fig. 22. — Estimated -Umin distribution for the typical timescale of t-^ = 40 days with (sohd 
hne) and without (dashed hne) the < 200 days cut. The same distribution with < 200 
days cut for the bright source star (/ < 17) events (dot-dashed hne), which is scaled to 
match the others at Wmin = 1, is also shown for comparison. 
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Mean Optical Depth <t{N^^^)> 

Fig. 23. — Probability P{t{N) > Tobs) that the optical depth t{N) in each simulated exper- 
iment is larger than the observed one Tobs as a function of the mean optical depth (T(A'exp)) 
of experiments in which an expected event number is A^exp- Ic confidence limits (dashed 
line) and the observed optical depth Tobs (solid line) are also presented. 
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Fig. 24. — Relation between tEm and the mean value of tEout for various source I-band 
magnitude as indicated in the figure. The mean < tEout > are taken over the events between 
each indicated / magnitude and / + 1 mag. We have no detection in ^Ein < 1 day for / = 22. 
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Fig. 25. — Example of the expected "OUTPUT" tEout distribution for tEin — 40 days (vertical 
dashed line), i.e., L'out(40). The mean < tgout > is 42 days in this case. 
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Mean INPUT Optical Depth <r^^(N^^p)> 

Fig. 26.— Relation between Mean optical depth of OUTPUT < rout(-/Vexp) > and INPUT 
< Tin(iVexp) > (solid line). Dashed line represents < rout(-^exp) >=< Tin{Ncxp) >■ The 
horizontal and vertical dotted line indicate observed (robs — 2.63 x 10~^) and estimated 
(r = 2.59 X 10~^) optical depths respectively. 
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Mean INPUT Optical Depth <Tin(^exp)> 

Fig. 27.— Probability P{Tout{N) > Tobs) that the OUTPUT optical depth Tout{N) in each 
simulated experiment is larger than the observed one Tobs as a function of the mean "INPUT" 
optical depth (rin(iVexp)) of experiments in which an expected event number is A^exp- 
confidence hmits (dashed line) and the estimated real optical depth r (solid line) are also 
presented. 
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Fig. 28. — Histograms of the event timescale tE distribution for 28 observed events with 
expected timescale distributions, normahzed to the observed number of events, for a fixed 
bar and disk density model and for various mass functions. The solid line indicates (i) Scalo 
(1986) 's PDMF; the dashed and dot-dashed lines indicates S function (ii) at M = O.IMq and 

(iii) at M = l.OM© respectively; the dotted and three dot-dashed lines show the power-law 

(iv) with a = 2.3, Mi = 0.1 and (v) with a — 2.0, Mi — 0.01, where the later one represents 
the brown dwarf rich mass function. 
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Fig. 29. — Histograms of the observed distribution, corrected for the detection efficiency. 
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Fig. 30. — Distribution of the contribution to the total optical depth (tE/e(^E)) of the 
observed event The contribution of the short timescale events to the optical depth is quite 
small. 
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Table 1. Selection criteria in cutl. 



type criteria 
type 1 Ndet,n = 

Ndet,p > 2 
^det,p 7^ ^clus,p 

Ratioi < 0.1 or Ratios < 0.1 
type 2 Ndet,p = 

■^det,n ^ ^ frame / ^ 
^clus,n ^ 5 

type 3 A^det,n 7^ and Ndet,p 

^clus,p ^ 4 
A^da6-,rt < 5 

{Ndet,p > 2 and Ndet,p ^ Nclus,p) or A^det,n > N frame/ 4: 



Table 2. Selection criteria in cut2. 
criteria 

Ntotal > 70 
Nin > 3 

Nout > 9 

1 < Np,ak < 3 
Ei,peafe > 20 

Nhi < 6 
xL < 4 or xln/xlut > 15 
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Table 3. Selection criteria in cut3. 



criteria 



X^/d.o.f. < 3.5 for Fpeak < 450, 000 ADU 
X^/d.o.f. < 100 for Fpeak > 450, 000 ADU 

^min ^ 1 

0.3 < < 200 day 



Table 4. Position of 28 microlensing event candidates 



field 


chip 


ID 


ID (real-time) 


ID (alert) 


RA (J2000) 


Dec (J2000) 


ngbl 


2 


2745 


- 


- 


17:58:13.136 


-29:09:14.23 


ngbl 


2 


4925 


2717 


2000-BLG-ll 


17:57:07.907 


-29:09:59.28 


ngbl 


2 


5076 


- 


- 


17:56:33.952 


-29:27:16.08 


ngbl 


2 


5157 


2667 


2000-BLG-7 


17:54:56.681 


-29:31:47.50 


ngbl 


3 


2567 


727 


2000-BLG-3 


17:54:29.770 


-28:55:59.31 


ngbl 


3 


6328 


2540 


- 


17:58:20.936 


-28:47:48.76 


ngbl 


3 


6344 


2548 


- 


17:55:05.425 


-28:50:34.60 


ngbl 


3 


7416 


- 


- 


17:57:54.728 


-28:54:32.66 


ngb2 


2 


3867 


1648 


- 


18:00:12.361 


-29:37:23.93 


ngb2 


3 


1932 


- 


- 


17:59:00.087 


-29:33:01.11 


ngb2 


3 


3807 


- 


- 


18:00:07.092 


-29:23:27.47 


ngb3 


2 


3465 


1316 


2000-BLG-9 


18:05:09.533 


-30:36:06.77 


ngb3 


3 


1041 


- 


- 


18:06:47.640 


-29:50:09.57 


ngb4 


1 


6678 


2806 


2000-BLG-13 


17:55:33.202 


-28:10:17.09 


ngb4 


3 


1293 


159 


- 


17:57:47.197 


-27:33:52.80 


ngbS 


1 


4316 


1673 


- 


18:01:06.711 


-28:52:22.28 


ngbS 


1 


4317 


1672 


- 


18:01:26.814 


-28:52:34.66 


ngbS 


1 


4318 


1668 


- 


18:01:44.791 


-28:58:03.53 


ngbS 


3 


1392 






18:01:21.393 


-28:02:51.70 


ngb6 


3 


2746 






18:04:46.138 


-28:31:31.54 


ngb6 


3 


3954 


1425 


2000-BLG-12 


18:03:54.776 


-28:34:58.62 


ngb7 


3 


2192 


703 


2000-BLG-8 


18:10:55.621 


-29:03:54.20 


ngb9 


3 


2336 


841 




18:10:17.990 


-27:31:19.31 


ngblO 


1 


1837 






18:08:32.453 


-26:09:29.66 


ngblO 


3 


2112 






18:08:51.926 


-25:24:40.46 


ngbll 


2 


1594 


1142 


2000-BLG-lO 


18:11:28.310 


-26:15:05.81 


ngbll 


3 


1063 






18:11:57.020 


-25:54:57.32 


ngbl2 


2 


3187 


1052 




18:14:47.421 


-25:32:53.64 



Note. — ID in this offline analysis, ID in real-time analysis and alert ID in Bond 
et al. (2001) are also presented. 
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Table 5. Parameters in microlensing light curve fitting for 28 candidates 



MOA ID 




fn 

''0 




A 






(Hi^Vl 




In 


field 


chip 


ID 


TD-245000 


1 owpr 

L\J VV V^l 






1 owpr 

L\J VV v^l 


npcif 






ngbl 


2 


2745 


1680.16 


200.0 


2.6 X 10^ 


8.8 X 10^ 


26.9 


41.3 


_ 


19.7 


ngbl 


2 


4925 


1799.38 


8.5 


8.6 


8.7 


56.9 


57.2 


57.5 


13.6 


ngbl 


2 


5076 


1725.17 


56.1 


91.7 


234.4 


33.2 


51.8 


106.4 


20.8 


ngbl 


2 


5157 


1725.79 


4.1 


5.1 


6.2 


2.9 


3.4 


4.0 


17.1 


ngbl 


3 


2567 


1691.82 


2.9 


9.4 


293.9 


6.9 


16.0 


_ 


17.5 


ngbl 


3 


6328 


1795.55 


9.4 


12.5 


16.8 


30.6 


37.2 


45.9 


17.6 


ngbl 


3 


6344 


1792.80 


14.9 


28.4 


60.3 


10.0 


15.4 


27.3 


18.3 


ngbl 


3 


7416 


1829.40 


10.6 


11.9 


13.2 


10.6 


11.4 


12.2 


16.2 


ngb2 


2 


3867 


1801.09 


4.6 


6.1 


9.1 


92.0 


116.2 


163.6 


17.4 


ngb2 


3 


1932 


1707.91 


76.5 


140.6 


361.4 


21.1 


36.3 


181.2 


19.9 


ngb2 


3 


3807 


1857.06 


7.2 


11.1 


17.6 


18.1 


25.6 


37.9 


17.8 


ngb3 


2 


3465 


1739.96 


28.7 


58.1 


115.9 


45.9 


60.3 


269.2 


18.8 


ngb3 


3 


1041 


1700.57 


1.2 


1.4 


1.9 


0.3 


7.0 


9.6 


15.1 


ngb4 


1 


6678 


1809.17 


16.4 


17.5 


18.6 


73.4 


77.8 


82.4 


16.2 


ngb4 


3 


1293 


1693.95 


3.0 


4.8 


12.8 


4.8 


5.5 


6.2 


14.5 


ngbS 


1 


4316 


1783.60 


33.5 


49.8 


89.0 


35.2 


45.6 


62.8 


18.1 


ngbS 


1 


4317 


1797.46 


3.8 


6.7 


16.9 


55.2 


85.2 


188.2 


17.9 




1 

X 


4318 


1778.51 


3.3 


5.0 


6.8 


22.0 


28.2 


34.2 


16.0 


ngb5 


3 


1392 


1672.15 


2.4 


4.1 


5.8 


1.6 


2.1 


3.0 


16.5 


ngb6 


3 


2746 


1714.69 


11.5 


16.5 


24.6 


16.3 


21.8 


30.9 


19.7 


ngb6 


3 


3954 


1796.00 


6.5 


7.5 


8.8 


12.4 


13.4 


14.5 


16.5 


ngb7 


3 


2192 


1732.84 


1.5 


1.7 


1.8 


12.7 


13.5 


14.4 


14.6 


ngb9 


3 


2336 


1792.14 


1.2 


1.4 


1.9 


0.8 


7.5 


10.3 


14.1 


ngblO 


1 


1837 


1663.93 


38.2 


51.8 


71.5 


12.7 


15.4 


19.0 


19.0 


ngblO 


3 


2112 


1700.85 


15.0 


28.8 




2.9 


5.8 


99.4 


18.3 


ngbll 


2 


1594 


1730.46 


9.0 


9.6 


10.4 


45.7 


48.8 


52.3 


16.6 


ngbll 


3 


1063 


1685.29 


16.5 


27.2 


82.4 


25.6 


38.4 


107.1 


19.7 


ngbl2 


2 


3187 


1789.45 


3.7 


4.6 


5.6 


6.7 


7.5 


8.2 


15.8 



Note. — J-band baseline magnitude of source star Jq is de-reddened to match the HST field 
by using the /-band extinction Ai map of each field (see §5). The symbol "-" means that a 
value could not be constrained. 



